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Abstract 



We show how parametric resonance mechanism follows quench mechanism in the 
classical linear sigma model. The parametric resonance amplifies long wavelength 
modes of the pion for more than 10 fm/c. The shifting from the quench mech- 
anism to the parametric resonance mechanism is described by a time dependent 
quantity. After the quench mechanism is over, that quantity has an oscillating 
part, which causes the parametric resonance. Since its frequency is 2m 7T (m n 
: pion mass ), very long wavelength modes such as k = 40MeV of the pion are 
amplified by the parametric resonance. 
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1 Introduction 



The notion of chiral symmetry breaking has a significant meaning to comprehend 
hadron physics. The lattice numerical simulations or theoretical studies predict that 
the chiral symmetry is restored under the extreme temperature. One of the purposes of 
experiments at Relativistic Heavy Ion Collider (RHIC) is to produce the chiral symmet- 
ric phase. With regard to those experiments, the idea of chirally misaligned vacuum, 
the disoriented chiral condensate (DCC) [0, 0, |3|, |, ^ , has been discussed. 

Rajagopal and Wilczek carried out numerical simulation in the classical 0(4) linear 
sigma model with fields = (a, 7?), which is often used to study the dynamics of 
chiral condensate ||. They show that, after quenching, long wavelength modes of the 
pion fields are amplified. This strongly suggests the formation of large DCC domains. 
Indeed, Asakawa, Huang and Wang verified later by numerical simulation that the 
DCC domains with 4 ~ 5 fm in size can be formed through a quench 0. If such 
large domains are actually produced in experiments, the existence of DCC would be 
confirmed by measuring, for example, a pion emission ratio R = n^o / (n^o + n v + + n n -) 

II- 

Therefore, it is crucial for the formation of large DCC domains that long wave- 
length modes of the pion are amplified after quenching. According to Ref. ||, one can 
understand the enhancement of such modes by considering the following mechanism. 
After quenching, the fields <fi roll down randomly in isospin space from the top of the 
Mexican hat potential. During the rolling down, the long wavelength modes of the pion 
field have negative mass square < 0, so that these modes grow exponentially. This 
is the mechanism Rajagopal and Wilczek argued, which we call quench mechanism. 

However, the numerical simulation in Ref. M shows that long wavelength modes of 
the pion are amplified continuously for more than 10 fm/c even when the effective mass 
becomes > after the rolling down. How can we understand such a phenomenon? 
It is necessary to answer this question since the formation of DCC domain is due to 
the enhancement of those modes. 

There are some attempts to explain such a phenomenon, which say that the oscil- 
lation of the sigma field around the true vacuum causes the enhancement of the pion 
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modes through parametric resonance [EL BL ITDJ. In these proposals, it is assumed that 



the field oscillates along the a direction around the minimum of its potential after 
the rolling down. It is, however, not obvious that the direction of the 0's oscillation is 
a, so that we do not adopt the assumption here. 

In this letter, we make a investigation of the cause of the enhancement of pion's 
long wavelength modes after the isotropic rolling down. The classical 0(4) linear sigma 
model is also used in our study, where we do not suppose that the field <fi oscillates along 
the a direction after the isotropic rolling down, as mentioned above. The point is that 
how the long wavelength modes of the pion are continuously amplified following the 
quench mechanism. In order to argue the quench mechanism, Rajagopal and Wilczek 
have analyzed the Euler-Lagrange equation for the pion modes with the approximation 
where the field <fi 2 (x) is replaced by its spatial average (0 2 )(t). (In Ref. |J, the spatial 
average (<j) 2 )(t) has been calculated by numerical simulation.) We also use the Euler- 
Lagrange equation with the same approximation to find how the pion modes develop 
after the isotropic rolling down. In our approach, the issue is the time dependence of 
the 0's spatial average (0 2 )(t), not the sigma field's one (o)(t) [@, pj, [ID]. The time 
evolution of the pion modes thus derived will be compared with the result of numerical 
simulation obtained by Rajagopal and Wilczek. 

In the next section, after reviewing the quench mechanism argued by Rajagopal and 
Wilczek ||, we describe how the parametric resonance mechanism follows the quench 
mechanism. The time evolution of the zero mode of the pion will be discussed in Section 
3. The last section is devoted to summary and discussion. 



2 Parametric resonance after isotropic rolling down 
of fields 

The linear sigma model is 

C = \d^d,<P - \\{<P 2 - v 2 ) 2 + Ha, (1) 

where = (a ,-n). Since an explicit chiral symmetry breaking term Ha exists, the 
pion has mass m 2 , = Hj and the true vacuum is = (f w ,0), where /„- satisfies the 
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following relation, 

- J) - H = . (2) 

Here the parameters are chosen as A = 20.0, v = 87.4 MeV and H = (119MeV) 3 so 
that 

U = 92.5 MeV, m w = 135 MeV, m a = 600 MeV. (3) 

In order to study the DCC, Rajagopal and Wilczek considered quenching in the linear 
sigma model. The quenching generates an isotropic rolling down of the fields <fi in 
isospin space from the top of the Mexican hat potential. They solved the equation of 
motion by a numerical simulation and found that long wavelength modes of the pion 
fields are amplified after quenching. The amplification is necessary for the formation 
of large DCC domains. To understand the enhancement of the long wavelength modes 
founded by the numerical simulation, they consider the equation of motion for the pion 
modes 7r(k, t) with the approximation where the field <p 2 (x) is replaced by its spatial 
average ((f) 2 ) (t) : 

^7?(k,t)+m e 2 ff (M)7f(M) = 0. (4) 

The effective mass m e s(k,t) depending on the time t and the momentum k = |k| is 
defined as 

mUKt) = k 2 + \{(<f)(t)-v 2 }. (5) 

The approximation, <fi 2 (x) — > (<f) 2 )(t), makes the equation of motion a linear equation 
with respect to the pion mode 7f(k, t) provided that (<f) 2 )(t) is given at first. During the 
isotropic rolling down of 0, the spatial average (0 2 ) is less than v 2 , so that the effective 
mass satisfies m 2 ff < for small k 2 . The long wavelength pion modes k with m 2 ff < 
then grow exponentially. This is the quench mechanism. 

What occurs after the quench mechanism is over? The effective mass becomes 
positive m 2 E > after the isotropic rolling down of <fi. We shall explain it in detail. 
Rajagopal and Wilczek calculated the spatial average (0 2 )(t) using the numerical so- 
lution of the (exact) equations of motion. According to Fig.2 in Ref. 0, (0 2 )(t) 
oscillates around f 2 (> v 2 ) with a frequency ~ 270 MeV = 2m n after the isotropic 
rolling down. The spatial average is then expressed approximately as 

(0 2 > (t) « f 2 + B 2 cos (2 m n t + 2 if), (6) 
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where the amplitude B 2 decreases with time. After B 2 < f 2 — v 2 , the spatial average 
(<p 2 )(t) is always larger than v 2 , so that the effective mass becomes m 2 s > 0. Therefore, 
it seems that all the pion modes are stable and oscillate after the isotropic rolling down 
(<f) 2 ) > v 2 , as discussed in Ref.||. However, the numerical simulation Fig.l in Ref. 
shows that the long wavelength pion modes are still amplified even when (0 2 ) > v 2 . 
Indeed, the pion mode with k = 0.20/a = 40 MeV is amplified until t ~ 28a ~ 28fm/c, 
where the parameter a = (200MeV) _1 . How can we understand this phenomenon? 

To understand the behavior of the long wavelength pion modes after the isotropic 
rolling down, we use the equation of motion (|]) with the same approximation, 2 (x) — > 
(<fr 2 )(t). As seen previously, the time dependent mass parameter m 2 ff (fc,t) in Eq.(f|) 
oscillates. When the oscillating part of m 2 ff (fc,i) has a suitable frequency, parametric 
resonance will work. Substituting Eq. @ into Eq. (f|), one has 

d 2 

— 7?(k, t) + {(ml + k 2 ) + \B 2 cos (2 m^t + 2 if)} vf(k, t) = 0, (7) 

where the term A B 2 cos (2m 7T t + 2ip) is called "stimulating force" ||. Before examining 
this equation in detail, let us consider it qualitatively. The eigenfrequency u n of the long 
wavelength modes of the pion is about , uj % = m 2 + k 2 ~ m n , and the frequency 
u s of the stimulating force is 2 , u s = 2 m n . Then, the parametric resonance works for 
the long wavelength modes, because u s m 2u} w [|TTJ. Consequently, the long wavelength 
modes of the pion with k m n are amplified by the parametric resonance after the 
quench mechanism is over. 

Here, our equation (0) should be compared with that of Mrowczyhski and Miiller 
|| to clarify the difference between them. They assume that, after the rolling down, 
the field <fi oscillates only along the sigma direction. Since the parametric resonance 
induced by the a field oscillation is concentrated, the equation for the pion modes 

7tSmOM) nas tne form B 1) 
d 2 

— 7fSL(k, t) + {(ml + k 2 ) - 2Ava cos (m a t + tp) } vf^(k, t) = 0, (8) 

where the frequency uj' s of the stimulating force is m a , u' s = m a . In this theory 

if the sigma's amplitude o"o is not large, the long wavelength modes of the pion with 

k <^ m w are not amplified by the parametric resonance because u' s ^> 2 u n . On the 
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other hand, in our theory (|?J), the long wavelength modes k <C m n are amplified because 
of uj s ~ 2 u n , as seen above. 

Let us go back to the discussion of our theory. To analyze Eq. ([?[) more quantita- 
tively, we introduce a new variable, z = m n t + (p . The equation (0) is then rewritten 

as 

d 2 

— fr (k, z) + {A + 2q cos (2z) } vf (k, z) = 0, (9) 

where 

A ml + k 2 \B 2 . . 

^ = 9=o-a. ( 10 ) 

2m£ 

and 5 2 decreases with time. For simplicity, we treat Eq. ([]) at short intervals of time, in 
which B 2 can be considered approximately as a constant. Eq.(||) is called the Mathieu 
equation || pT] , whose unstable solution increases rapidly with time. The unstable 
solutions 7? (k, z) of the Mathieu equation in the first resonance satisfy the following 
condition \A — 1| < q + 0(q 2 ) if < A < 2 and q < 1. Therefore, the long wavelength 



modes with 



fc<-Bi/^, (11) 



are amplified by the parametric resonance. 



Does the condition (|TTJ) explain quantitatively the enhancement of the long wave- 
length pion modes observed in the numerical simulation of Ref. || ? We shall focus 
on two stages of the parametric resonance: (i) at the early stage and (ii) at the late 
stage, (i) At the early stage with the amplitude B 2 = f 2 — v 2 = ml/X ~ (30MeV) 2 
corresponding to the time t ~ 10a ~ lOfm/c in Fig. 2 of Ref.||, the space average 
{4> 2 )(t) oscillates as v 2 < (<p 2 )(t) < 2f 2 — v 2 . The condition ( |TTD predicts that the pion 
modes with 

Tn 

k < —jL = 96 MeV, (12) 
v 2 

are amplified when t ~ lOfm/c. This prediction is verified if one compare it with 
the result of the numerical simulation Fig.l of Ref.0. (ii) At the late stage with the 
amplitude B 2 ~ (15MeV) 2 corresponding to the time t ~ 26fm/c in Fig. 2 of Ref.||, 
(4> 2 ){t) oscillates as v 2 < (4> 2 )(t) < 2f 2 - v 2 . This time the condition (TTTj) predicts 

k < 47 MeV . (13) 
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This prediction is also verified in the same manner. Thus the parametric resonance (|TTj) 
is able to explain quantitatively the enhancement of the long wavelength modes of the 
pion < k < m n observed in the numerical simulation ||. The characteristic features 
of the enhancement after the quench mechanism are (a) that the long wavelength modes 
are amplified for more than 10 fm/c and (b) that the lower the momentum k of the 
pion is, the longer a span of the enhancement becomes. The time evolution of the zero 
mode k = of the pion will be discussed in the next section. 

3 Pion's zero mode oscillation 

In the parametric resonance mechanism, the oscillating part of (4> 2 )(t) plays a role of 
the " stimulating force" . What is the origin of the oscillation part B 2 cos (2 m n t + 2 if) 
having the frequency 2 m„ ? The spatial average (<f) 2 )(t) is written by the sigma modes 
cr(k, t) and the pion modes 7r a (k, t) as 

^ = I^m\ 2 +±J^\m 2 - < 14 > 

According to the numerical simulation in Fig.l of Ref.||, the frequency of |7r a (k, t)\ 2 
is given by 2^Jm 2 + k 2 , while the frequency of |a(k, t)| 2 can not be specified due to its 
complicated behavior. Therefore, the (</> 2 )'s oscillating part B 2 cos (2 m n t + 2 if) arises 
from the zero mode k = of the pion. 

Now, we suppose that the oscillation of the pion's zero mode is represented as the 
following simple form 

(7r )(t) = 7 a cos(m 7r t + <p), (a = 1,2,3), (15) 

in order to find the relation between B 2 and the amplitude 7 a of the pion's zero mode 
oscillation. We obtain 

1 3 

J B 2 cos(2m^ + 2<^) = - ^ ll cos ( 2m nt + 2<p), (16) 

because the origin of the oscillating part of {<p 2 )(t) is the zero mode of the pion. Con- 
sequently, it yields a relation 

B 2 = \Y.ll- (17) 
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For example, if the zero mode has a form (v?)(£) = (7, 7, 7/4) cos (m n t + ip), the relation 



fllT]) becomes B = y 33/32 7. Let us estimate the pion's zero mode amplitude 7 at the 
two stages of the parametric resonance (i) and (ii) discussed in the previous section, 
(i) At the early stage ( t ~ lOfm/c, B 2 = f 2 — v 2 ), the magnitude of the amplitude is 
7 ~ 0.32 U (h) At the late stage ( t ~ 26fm/c, B 2 = (15MeV) 2 ), that is 7 ~ 0.16 f n . 

According to the numerical simulation Fig. 2 of Ref. 0, the amplitude B 2 decreases 
with time. The relation fllTD implies that the amplitude 7 a of the pion's zero mode 
also decreases. This leads us to the following problem. The fact 7 a decreases with 
time contradicts with the statement in the previous section: the statement is that the 
zero mode k — of the pion field belongs to the unstable solutions ( |TTD of Eq.(^), 
so that the zero mode is amplified. This contradiction is, however, resolved since the 
approximation used in Eq.(^) is not applicable to the zero mode. The reason is as 
follows. The equation of motion ([^D with the approximation is regarded as a linear 
equation with respect to the pion mode ir a (\s.,t). Hence the nonlinear terms of 7r a (k, t) 
involved in the exact equation of motion are not concerned. But, in the equation of 
motion for the zero mode k = 0, the effect of nonlinear terms of zero mode is not 



negligible JLO 



4 Summary and discussion 

We studied the classical linear sigma model to investigate the disoriented chiral conden- 
sate ( DCC ). We showed that the long wavelength modes of the pion are amplified for 
more than 10 fm/c by the parametric resonance. The parametric resonance mechanism 
follows the quench mechanism, both of which enhance the long wavelength modes of the 
pion. The shifting from the quench mechanism to the parametric resonance mechanism 
is described by a time dependent quantity. That is the spatial average ((fi 2 )(t): 



{<f){t) 



< v 2 , early time 

f 2 + B 2 cos (2 m n t + 2 ip) , late time . 



During (</> 2 )(t) < v 2 , the quench mechanism works as argued by Rajagopal and Wilczek. 
After that, (</> 2 )(t) oscillates around f 2 with the frequency 2m n . This oscillation plays 
a role of "stimulating force" for the pion modes k(> 0). Since the frequency of the 
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stimulating force is about twice as much as the eigenfrequency of the long wavelength 
pion modes, the parametric resonance mechanism works. That mechanism can explain 
quantitatively the result of numerical simulation in Ref. concerning the amplification 
of the long wavelength modes of the pion. 

As stated in the previous section, the "stimulating force" with the frequency 2m 7T 
is attributable to the zero mode of the pion. Therefore, we believe that the oscillation 
of the pion's zero mode k = around the true vacuum causes the parametric resonance 
by which the long wavelength modes k (> 0) of the pion are amplified. This scenario is 
contrast with the scenario of Refs. || |9], [T(J in which the oscillation of the sigma's zero 
mode causes the parametric resonance. 
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